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The  research  reported  upon  here  addresses  several  elements  of  the  structure  of  nearly 
parallel  turbulent  shear  flow.  Both  boundary  and  free  shear  layer  flows  In  constant  as 
well  as  variable  pressure  fields  are  addressed,  with  respect  to  variable  pressure  condi¬ 
tions,  unsteady  as  well  as  steady  flows  are  embraced  by  the  approaches.  As  a  serendipi¬ 
tous  consequence,  many  of  the  results  also  apply  to  flowfield  separation  in  time-dependent 
flow,  a  topic  of  current  interest  to  AFOSR.  Specifically,  new  (and  soon  to  be  published) 
aspects  of: 


•  Turbulent  structure  advection  velocities, 

•  Turbulent  boundary  layer  "instability",  and 

•  Turbulent  shear  flow  angular  structure 


have  been  exposed.  Due  to  the  inherent  "try  and  fail"  nature  of  basic  research,  the 
approximately  one  man-year  effort  invested  so  far  has  not  concluded  the  study  of  any 
of  these  three  areas,  although  four  publications  have  resulted  from  the  investigation. 
The  work  has,  however,  established  the  feasibiiity  of  further  inquiry  benefiting: 


•  Basic  turbulence  physics, 

•  LES  and  turbulence  modeling,  and 

a  Flowfield  separation  physics  In  steady  and  unsteady  flow. 

In  particular,  a  coordinated  experimental  and  computational  study  of  those  turbulent 
shear  flow  details  which  have  been  exposed  would,  in  fact,  serve  the  dual  goals  of 
understanding  turbulence  structure  as  well  as  separation  physics. 
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As  a  direct  result  of  the  work  reported  here,  an  SBIR  study  of  turbulence  in  hypersonic 
flow  was  initiated  for  AFOSR,  the  fruits  of  which  now  support  an  advanced  R&D  effort 


at  Eglin  AFB. 
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Discussion 


The  focus  of  the  research  reported  upon  here  is  the  structure  of  nearly  parallel  turbu¬ 
lent  shear  flows  in  constant  and  variable  pressure  fields.  Both  wall-bounded  and  free 
shear  layers  were  studied,  as  appropriate,  and  the  variable  pressure  results  apply  to 
unsteady  as  well  as  steady  flows.  In  addition,  many  of  the  results  bear  directly  upon  the 
generalized  separation  problem,  addressed  here  insofar  as  it  is  relevant  to  turbulent 
structure  details.  This  is  reminiscent  of  Prandtl's  view  that  turbulence  is  sustained  by 
local  separations. 

The  work  discussed  here  is  scheduled  to  be  submitted  for  publication  during  the  next 
few  months.  In  addition,  various  other  tangential  papers  which  have  resulted  from  this 
work  have  been  published  and  are  included  in  Appendix  A. 

In  reporting  upon  the  achievements  of  this  research,  it  is  easy  to  present  them  lineariy 
in  an  orderly  and  coherent  form.  In  actuality,  the  procedures  Involved  where  chaotic  and 
not  serial  but  strongly  linked  parallel  efforts  that  took  form  only  after  considerable  stum¬ 
bling  and  "probing  for  soft  spots"  in  the  mathematics.  Therefore,  with  this  admission  in 
the  forefront,  the  discussion  begins  arbitrarily  with  a  report  on  work  carried  out  relative 
to  the  propagation  of  large  scale  structures  in  turbulent  shear  flows.  The  only  rationale 
for  beginning  with  this  is  its  simplicity  and  possible  surprisingly  large  future  impact. 

I.  Advectlon  Velocities 

The  following  generalized  evolution  equation: 


arises  in  several  applications  relevant  to  this  work: 

1.  As  an  integral  boundary  layer  thickness  equation, 

2.  As  the  equation  for  an  instability  envelope  in  a  two-time 
perturbation  analysis, 

3.  As  the  governing  equation  for  the  flow  angle,  and 

4.  As  a  pdf  transport  equation. 

The  second  and  third  source  of  Eq.  (1)  shall  be  revisited  in  connection  with  the  next 
aspects  of  the  work  to  be  reported.  No  advantage  has  yet  been  taken  of  the  congruence 
to  the  third  area. 

For  our  purposes,  the  (two-dinr^ensional  or  spanwise-averaged  three-dimensional)  time- 
dependent  equations  of  motion  may  be  integrated  across  the  flow*  to  provide  Eq.  (1), 
wherein  the  following  definitions  apply: 

L  =  0 ,  the  momentum  thickness, 

U^=  U/H,  the  freestream  velocity  divided  by  the  shape  factor, 

P  =  generalized  streamwise  gradients  term,  and 

G  =  C*/2  ♦  (v  AJXI  -  the  skin  friction  coefficient 

plus  vectored  wall  toss  transfer  terms,  times  U  . 

Of  interest,  initially,  was  simply  how  fast  large  scale  structures  propagate.  Clearly, 
from  Eq.  (1),  a  characteristic  speed  equal  to  H"*  times  the  local  freestream  velocity 
appears.  For  constant  pressure  turbulent  boundary  layers,  H  «  1.4  and,  therefore,  H**« 
0.7,  consistent  with  experimental  evidence  re^rdlng  certain  of  the  largest  scale  struc¬ 
tures.  For  our  purposes,  therefore,  H  Is  now  defined  by  the  expr^on: 
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f 


(2) 


H  =  UAJ^ 
c 


whereas,  it  is  computed  from  6*/0,  the  ratio  of  displacement  to  momentum  thicknesses, 
with  respect  to  interpreting  the  equations,  their  solutions  and  H,  physically,  this  orienta¬ 
tion  is  more  useful  than  discussing  the  6*/e  scale  ratio. 

In  order  to  solve  Eq.  (1)  for  the  evolution  of  a  specific  burst-like  structure  in  the  shear 
flow,  P  must  be  computed,  ostensibly  from  the  terms  in  its  defining  equation.  For  this 
study,  however,  an  inverse  approach  was  taken  in  that  P  itself  was  specified  parame¬ 
trically.  Subsequently,  the  definition  for  P  was  used  to  solve  for  the  implied  and  H 
distribution  from  the  equation: 

(|t  »  =  P  .  (3) 

Of  prime  importance  here  is  the  appearance  of  two  additional  characteristic  speeds  of 
propagation, 


u^'/u  =  (H  ♦  2)/2H  ,  and 

u  ••  =  2U^’ 
c  c 


(6) 


For  a  constant  pressure  turbulent  boundary  layer  with  H  *  1.4,  both  u  '  and  U  *•  are 
considerably  greater  than  the  freestream  velocity,  consistent  with  experimental  obser¬ 
vations.  In  addition,  several  aspects  of  these  functional  forms  for  the  three  characteris¬ 
tic  speeds  are  clarified  by  the  following  plot: 
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Over  the  range  of  data  encountered  experimentally,  both  u  and  U  ”  are  always  less 
than  and  greater  than  freestream,  respectively,  for  both  laminar  and  turbulent  flows. 

For  laminar  flow,  u^*  is  always  less  than  freestream.  However,  for  turbulent  flow  u_' 
transitions  from  greater  than  freestream  to  less  then  freestream  velocity  at  Hs2.  It  is 
the  hypothesis  of  this  work,  therefore,  that  a  significant  structural  change  in  the  flow 
relative  to  the  propagation  of  large  scale  structures  occurs  at  this  point.  Ciearly,  beiow 
a  value  of  H=2.0,  the  freestream  is  not  able  to  feed-forward  information  on  H  (carried 
by  u^'),  as  can  occur  above  Hs2.0.  In  addition,  the  broad  band  of  experimental  H  values 
that  characterize  turbulent  separation  (~1.8  *  2.3)  may  merely  reflect  wave  dynamics 
peculiar  to  the  particular  pressure  distributions  whereas  the  critical  value  of  H  is  indeed 
2.0. 
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with  the  aid  of  semi-empirical  H  and  laws*,  a  P  function  which  linearly  ramped  up 
in  time  to  a  Gaussian  spatial  distribution  was  input  to  a  numerical  solution  of  Eqs.  (1) 
and  (3).  Two  peak  magnitudes  of  P  were  identified  which  forced  the  peak  H  to  be  in 
the  range  of  2.0  and  somewhat  smaller,  respectively.  When  the  flowfield  solution  attain¬ 
ed  equilibrium,  a  disturbance  was  injected  into  the  flow  in  the  form  of  pulsed  and  vec¬ 
tored  mass  transfer.  For  the  smaller  of  the  two  peak  P  and,  therefore,  H  levels,  the 
disturbance  was  not  catastrophic.  However,  at  the  larger  peak  P,  the  disturbance  caused 
the  shape  factor  to  exceed  2.0  and  grow  exponentially,  as  shown  in  the  accompanying 
Figs.  1&2,  a-c.  These  results  also  point  out  the  effect  of  this  inverse  approach  wherein 
the  local  external  flow  responds  to  the  presence  of  the  'bulge'  in  the  boundary  layer. 
This  effect  is  critical  to  exposing  the  importance  of  H=2.0. 

In  conclusion,  to  the  extent  possible  with  the  time  and  funds  available,  the  importance 
of  U  *  and  its  behavior  as  H  changes  in  space  and  time  due  to  steady  and  non-steady 
pressure  gradients  has  been  established. 


Figs.  1  &  2,  a-c. 

Figs.  1:  Damped  Response  of  TBL  to  Pulsed  Blowing.  H, 

and  v^  =  0.003.  undisturbed 

la.  Momentum  Thickness 

lb.  Freestream  Velocity 

lc.  Shape  Factor 

Figs.  2;  Undamped  Response  of  TBL  to  Pulsed  Blowing.  K  i'97 

and  mJaj  .  0.003.  undisturbed 

w 

2a.  Momentum  Thickness 
2b.  Freestream  Velocity 
2c.  Shape  Factor 
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Figure  la.  CONTOUH 
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Figure  Ic.  jSSk 
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Figure  2a. 
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Figure  2c.  A 
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II.  Flow  Angle  Structure 


A.  Background 


Previous  discussion  relative  to  pulsed  and  vectored  wall  mass  transfer  indicated  how 
this  research  first  stumbled  onto  the  importance  of  flow  angles  relative  to  vorticity 
transport.  In  particular,  the  vectored  wall  mass  transfer  term: 


(1  -  u^/U)  (v^/u) 


used  to  generate  the  pulsed  disturbance  discussed  in  the  part  I.  of  this  section,  was  stu¬ 
died  in  the  first  paper  included  in  Appendix  A.  Written  in  terms  of  the  variables: 


«  =  and 


the  vectored  mass  transfer  term  becomes: 


{m*-m(l+o* )^}  {-o/(l+o*)} 


Equation  (7)  exhibits  the  following  extrema: 


o/U  *  1/2,  for  a  given  o,  and  • 


u  /u  *  (1±(1  ♦  8m*)*)/4,  for  a  given  m. 


15 


Subsequently,  the  research  led  to  the  second  paper  in  Appendix  A,  wherein  the  detaiis 
of  vectored  and  pulsed  mass  transfer  on  the  boundary  conditions  are  presented.  The 
essence  of  that  work,  as  it  relates  to  the  angular  structure  of  turbulent  shear  flows, 
may  be  presented  in  terms  of  the  two-dimensional  incompressible  Navier-Stokes  equations: 

H  +  VO)  and 

X  y 

(9) 

H  -  voj  , 
y  X  ’ 

where:  0*=  v/u  , 

cj  s  vorticity  and 

H  =  p/p  +  (u*+v*)/2 

Note  that,  evaluated  at  the  mass  transfer  surface,  these  equations  prove  that  the  vorti¬ 
city  flux  varies  even  in  constant  pressure  flows  if  the  mass  transfer  is  vectored,  pulsed 
or  exhibits  spatial  variation. 

By  using  the  condition  of  flowfield  incompressibility,  Eqs.  (9)  may  be  cross-differentiated 
to  provide  the  following  expression  for  H: 

7*h  =  u*  ♦  u  (u)y  -  00)^)  .  (10) 

Inasmuch  as  streamlines  (r|>)  are  defined  by:  dil^udy-vdx,  along  4*  «  constant  lines,  dy/dx  s 
6.  Therefore,  the  H  field  responds  to  gradients  in  vorticity  normal  to  the  streamlines. 

*Note  that  6  was  also  used  to  represent  the  momentum  thickness  in  part  I. 


u^+6(uaj)  =  - 
v^  -  1  (uoj)  =  - 


By  manipulating  Eqs.  (9),  the  following  expression  ensues: 

0.  -  CO  =  -  (H  -  0H  )/u  -  v(a)  +  0(0  )/u  ,  (11) 

L  y  ^  ^  y 

where  ( 1+0* )  has  been  approximated  by  1  in  the  coefficient  of  u)  for  weakly  non-paral¬ 
lel  flows. 

The  interpretation  of  Eq.  (11)  is  that  0^  is  a  measure  of,  but  lags  or  leads  the  vorticity 
according  to  cross-streamline  gradients  in  H  and  along-streamline  gradients  in  u. 

with  this  motivation,  the  basic  angular  structure  of  turbulent  shear  flows  is  being  analy¬ 
zed  to  expose  the  physics  and  aid  the  modeling.  As  this  report  is  being  written,  a  fully 
three-dimensional  solution  of  a  nominally  two-dimensional  free  shear  layer  is  being  com¬ 
pleted,  which  forms  the  basis  for  a  future  investigation  and  publication. 


B.  Development 


For  the  purposes  of  exposition,  consider  initially  a  two-dimensional  flowfield  ( 
and  examine  the  behavior  of  the  angular  tangent  variable : 

0  =  v/u 

By  applying  traditional  Reynolds  decomposition,  it  may  be  shown  that: 

^  -  0  =  -  8 • u ' /u  , 

where:  6  =  v/u 

In  addition: 

& ' u '  =  ( u ' v  *  -  Su ' * ) /u 

By  combining  Eqs.  13  &  14,  the  following  expression  ensues: 

S  -  ©  s  R  ♦  5(u' ^/u* )  » 

where:  R  =  u*v»/u*  , 

herein  referred  to  as  the  local  Reynolds  stress  correlation  coefficient. 


u, v;x,y) 


(12) 


(13) 


(14) 


(15) 


For  small  5  and  turbulence  intensity,  Eq.  (15)  nnay  be  approximated  by: 


(16) 


which  has  the  obvious  interpretation  that  the  difference  between  the  mean  flow  angle 
and  the  angle  of  the  mean  flow  is  equal  to  the  local  Reynolds  stress  correlation  coeffi¬ 
cient.  In  addition,  even  on  lines  of  symmetry  when  0  =  0,  S  may  be  non-zero  and  is 
equal  to  ^  with  the  aforementioned  restrictions. 

In  order  to  form  various  correlations,  neglecting  higher  order  fluctuation  products,  the 
fluctuating  angle  may  be  written: 


0'  =  ©(v'/v)  .  5(u’/u)  ♦  H.O.T. 


(17) 


where: 


H.O.T.  =  (STTJT  -  e'u')/u 


In  the  four  quadrant  decomposition’  of  coherent  structure  data,  the  second  and  fourth 
quadrant  contribute  most  to  the  net  production  of  Reynolds  stresses,  associated  with 
the  ejection  and  sweep.  Clearly  from  Eq.  (17),  e*  is  largest  when  v’  and  u*  are  of  the 
opposite  sign,  as  in  the  second  and  fourth  quadrant.  However,  v*  and  u*  are  weighted 
by  0  and  which  differ  by  R7  Therefore,  the  behavior  of  S  and  0'  is  significant  and,  in 
conjunction  with  conventional  approaches  to  obtain  v  and  u,  and  thus  6 ,  may  be  a  mwe 
effective  way  of  exposing,  e.g.,  the  Reynolds  stress  structure. 


f 
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C.  Analysis 


As  shown  in  Appendix  B,  the  full  N-S  equations  may  be  written  in  terms  of  angular 
variables  in  the  form: 


where: 


(18) 


Njj.  =  a/at  u  (a/ax  +  Ga/ay  +  a/az)  , 

0,a  =  v/u,  w/u  , 

F,G,H  =  -  +  vV*u.  ;  i  =  1,2,3  ,  and 

X.  1 

P  =  p/p. 

Consider,  for  simplicity,  an  essentiaiiy  inviscid,  two-dimensionai  flow  such  that  Eq.  (18) 
and  the  continuity  equation  become: 

®t  * 

(19) 

“x  *  ^“®^y  *  ® 


20. 


Written  in  streamf unction  variables,  Eq.  (19)  becomes: 


+  uG^  =  0  (P^  -  u0P^)/u  .  (20) 

By  expanding  0  in  a  series,  the  first  order  weakly  non-parallel  solution  obeys  the  follow¬ 
ing  equation: 


0  +  u0  =  0  P  /u  +  H.O.T. 

w  X  X 


(21) 


Conveniently,  this  is  the  precise  form  of  the  momentum  thickness  evolution  equation 
studied  in  part  I.  of  this  section.  We  have  also  exposed  a  new  parameter,  P^^/u,  to  inves¬ 
tigate  such  that  the  understahding  of,  and  ability  to  model  coherent  structure  is  re¬ 
inforced.  As  this  report  is  being  written,  a  fully  three-dimensional  N-S  solution  of  a 
nominally  two-dimensional  planar  mixing  layer,  such  as  that  shown  in  Figs.  3,  is  being 
completed  for  fluture  publication,  along  with  a  refined  version  of  Appendix  B.  The  fluc¬ 
tuation  data  obtained  shall  be  used  to  evaluate  the  angular  fluctuations  and  various 
correlations  thereof  with  respect  to  turbulent  burst/sweep  physics  and  transport  modeling. 
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The  pioneering  work  of  Landahl"  in  weakly  non-linear  instability  (wave-iike)  representa¬ 
tions  of  turbulence  forms  the  basis  for  the  computations  of  coherent  structure  fluctua¬ 
tions  which  have  been  carried  out.  A  multiple-elennent  decomposition*'^  of  the  governing 
equations  was  used  in  which  temporally  coherent  and  two  distinct  spatial  scales  of  tem¬ 
porally  incoherent  fluctuations  are  'superposed'  in  order  to  constitute  the  solution.  By 
judicious  use  of  appropriate  time,  phase  and  spatial  averaging,  governing  equations  for 
each  element  may  be  derived,  as  shown  in  Fig.  4. 


A  FOUR-ELEMENT  DECOMPOSITION  OF  THE  NAVIER-STOKES  EQUATIONS 


*"1  ♦  „  ‘"i .  1  ^  “'‘'1 

■sr  "j  ■jr  *  p  -sq"  *  *  arsT 


Decompose  the  flow  Into  four  components  by: 


^1  «  ♦  IT^  +  U]  +  uj 


p  •  P  *  P  *  P*  -fp'  , 


where: 


,  P  ■  time  mean  flow 

(Tf,  P'  •  coherent  component 

Uj,  P'  «  Incoherent  component.  Urge  scale 

u{»  P*  *  Incoherent  component,  smII  scale 
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where  -  ■  time  average 

<  >  •  phase  average 

{  }  *  large  scale  average 


Define  fluid  stresses: 

'■ij  *  “i“j  • 

Further,  decompose  ^ese  streses  as: 

*^1j  "  ^  *^ij  *  -1j 

*  ^3  *  • 
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Of  primary  interest  in  the  research  reported  here  is  the  large-scale  temporally  incohe¬ 
rent  fluctuating  element,  generally  referred  to  as  (spatially)  coherent  turbulent  structure. 
The  small-scale  temporally  incoherent  element  enters  into  and  affects  the  coherent 
structure  solution,  however.  Rather  than  solving  the  small-scale  equations,  the  terms  in 
the  large-scale  equations  resulting  therefrom  are  modeled  utilizing  classical  gradient 
transport  theory.^’®’® 

The  virtue  of  retaining  temporally  coherent  elements  in  the  decomposition  is  that  un¬ 
steady  flow  effects  on  the  turbulent  structure  may  be  readily  incorporated.®  In  parti¬ 
cular,  the  effect  of  various  time-dependent  freestream  velocity  fields  on  the  turbulent 
transport,  insofar  as  it  affects,  for  example,  flowfield  separation,  could  be  addressed 
by  this  formulation. 

The  conventionai  wisdom  regarding  the  constant  pressure  mean  turbulent  flow  velocity 
profile  is  that  it  is  stable.  This  is  clear  form  the  work  of  Landahl  in  which  fluctuation 
correlations  arising  from  non-linear  terms  are  treated  as  a  known  non-homogeneous  func¬ 
tion.  The  resultant  Orr-Sommerfeld  problem  generates  no  unstable  modes.  In  the  present 
work,  the  temporally  incoherent  small-scaie  processes  in  the  large-scale  equations  are 
considered  to  be  (at  least  partially)  functions  of  the  large-scale  solution.  The  modeling 
employed  is  gradient  transport  with  a  spatialiy-variable  and  non-equilibrium^®  (small- 
scale  eddy  viscosity)  coefficient,  c.  The  non-equilibrium  is  specified  to  be  of  the  form: 

De/Dt  ♦  €  s 
e  eq 


► 


(22) 


The  objective  of  tMs  aspect  of  the  work  was  to  review  the  stability  characteristics  of 
the  mean  turoulent  velocity  profile  with  a  spatially-variable  and  non-equiJibrium  small- 
scale  process  mode!  and  to  structure  the  solution  for  spatially-  and  temporaliy-variable 
pressure  fields.  The  transformed  weakly  non-linear  vorticity  equations  are  shown  in  Fig. 
5.  Note  also  the  possibility  of  resonance  between  the  vertical  vorticity  and  vertical 
velocity. 


Fig.  5. 


FOURIER  AND  LAPLACE  TRANSFORMED  3-D  VORTICITY  EQUATIONS 
LINEARIZED  LARGE  SCALE  MOTIONS 


|{U-C)  ♦  1  [(Re‘^+  c)(02-k2)  ♦  -|(Du)  ♦!  (Oc) (2D2-k2)|i, 

-  1“  |(Dc)0  -  (2D*«)|u2 

|{u^c)  ♦  1  [(R®*^+  e)(02-k*)  ♦  (Dc)D]|w2  ■  -  ^  (Dir)u2 


|{u^c)  ♦  1  [(R®*^+  e)(02-k*)  ♦  (Dc)D]|w2  ■  -  ^  (Dir)u2 

|(u-c)  ♦  i  [(Re'^  f)(D2-k2)  ♦  (02c)^i3.  |(oJ)  ♦!  (0e)(2D*-k*)|ilj 

♦  {(Dc)D  -  (2D2o|u2 

♦  {-  i  [0^  ♦  (0u)D]|m2 
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The  mean  velocity  fields  presented  in  Appendix  C  were  taken  from  NBS  data  due  to 
Klebanoff  and  a  wide  range  of  values  for: 

6  =  real  span-wise  wave  number, 
w  =  real  frequency, 

T^=  non-equilibrium  time  constant,  and 
=  peak  value  of  e/v  , 

mag 

were  investigated.  (Note  that  non-dimensionalization  with  respect  to  6*  and  U  is  assumed 
throughout).  The  eddy  viscosity  distributions  which  were  utilized  are  also  shown  in  Appen¬ 
dix  C. 

As  with  Landahl's  work,  variations  in  the  solution  with  B  and  w  were  relatively  un-dra- 
matic.  Therefore,  at  6=0.1  and  u)=2.0,  a  wide  range  of  values  was  studied.  Surpri¬ 
singly,  the  effect  of  non-equilibrium  on  the  solutions  was  relatively  minor.  At  a  of  1.0, 
various  values  and  c  spatial  distributions  were  also  studied.  As  long  as  a  spatially 
non-uniform  e  of  the  general  form  shown  previously  was  utilized,  unstable  wave  modes 
were  observed.  The  typical  results  presented  here  in  Fig.  6  indicate  that,  with  appropriate 
small-scale  process  modeling,  linear  instability  (indicated  by  Ksa^/a^  <  0)  of  even 
the  constant  pressure  mean  turbulent  flow  velocity  profile  is  possible.  Additional  results 
are  presented  in  Appendix  C. 


ALPHA  IMAGINARY 


The  Hokenson  Company 


3.00 


2.00 


1.00 


0.00 


Fig.  6. 
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Inclusion  of  non-paraliei  effects  anc  (spatially  and/or  temporally)  variable  freestream 
velocity  effects  are  included  utilizing  conventional  multiple-scale  perturoation  techni¬ 
ques.  Note  that  both  spatial  and  temporal  pressure  variations  are  represented  by  the 
temporally  coherent  element  of  the  decomposition  for  convenience.  In  addition,  these 
effects  and  this  procedure  is  identical  to  that  carried  out  for  the  evaluation  of  weakly 
non-linear  terms  in  the  large-scale  structure  solution  itself,  in  the  absence  of  temporally 
coherent  effects.  In  each  of  these  cases  an  amplitude  evolution  equation  which  envelopes 
the  oscillation  is  derived  and  found  to  be  identical  in  form  to  that  of  the  generalized 
shear  layer  bulqe  evolution  and  anqular  structure  equations  previously  studied. 


Result?  and  Conclusions 


As  a  result  of  this  research  carried  out  under  contract  to  AFOSR/NA,  the  following 
discoveries  may  be  enumerated: 


•  In  non-steady  turbulent  shear  flows,  one  of  the  dominant  advection  velocities 
transitions  from,  'super'  to  'sub*,  relative  to  the  local  freestream  velocity,  at 
a  value  of  the  shape  factor  (h)  equal  to  2.  At  this  point  the  freestream  is 
able  to  feed-forward  information  regarding  the  disturbance  and  catastrophic 
growth  leading  to  separation  occurs. 

•  The  velocity  vector  angle  structure  of  turbulent  shear  flows  appears  to  be  a 
useful  variable  for  interpreting  coherent  structure  data,  predicting  its  evolu¬ 
tion  and  formulating  turbulent  transport  models. 

•  If  adequate  models  of  the  apparent  stresses  in  constant  pressure  turbulent 
boundary  layers  are  employed,  linearly  unstable  modes  arise.  Extension  to 
spatially-  and/or  temporally-variable  pressure  fields  is  straightforward. 


Based  on  these  accomplishments,  it  is  concluded  that  further  work  is  warranted  relative  to: 


a  An  experimental  and  detailed  numerical  study  to  corroborate  the  observed 
criticality  of  H  r  2. 

a  Fully  three-dimensional  numerical  sinxilations  and  experimental  investiga¬ 
tions  of,  for  example,  free  shear  layers  in  terms  of  their  angular  structure, 
relative  to  both  coherent  structure  interpretation  and  modeling. 

a  The  solution  of  a  full  initial  value  problem  utilizing  the  weakly  non-linear 
wave-guide  simulations  of  large-scale  structure. 


The  final  result  and  conclusion  of  this  effort  is  that  it  motivated  an  SBIR  study  on 
hypersonic  turbulence  for  AFOSR,  the  fruits  of  which  now  support  and  advanced  R&D 
project  at  Eglin  AFB. 
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The  exceUent  data  presented  in  Ref.  I  provide  substan¬ 
tial  guidance  in  the  understanding  of  many  of  the 
subtleties  of  turbulent  boundary-layer  (TBL)  flows  with 
blowing,  and  the  modeling  thereof.  In  an  effort  to  peel  away 
quantitatively  the  various  interacting  elements  of  the  flow,  it 
is  proposed  that  the  following  approach  may  be  useful  for 
exploiting  such  data  in  model  devdopment.  If  the  Reynolds- 
averaged  equations  of  motion  are  further  locally  spatially 
averaged  in  horizontal  planes,  with  a  scale  large  rdative  to 
the  porosity  yet  small  r^tive  to  tlw  boundary-layer  growth, 
the  ditcrete  jet  btowing  may  be  converted  into  in  distribnied 
mast  uantfer  equivalent.  For  example. 


iiw<ii><fn*  (1) 

where  n  it  the  tfane-avenged  ttceaniwte  velocity  and  <  ), 
(  )*  denote  the  local  spatial  avenge  and  the  variatieo 
therefrom,  tespeedvdy.  As  a  lesiift  of  this  local  spatial 
averaging,  the  governing  aquation  ntmlincaritiet  generate  ex¬ 
tra  (Reynolds  stress-like)  terms  due  to  die  spatial  vartetkmt. 
(Note  that  a  piecewise  oonstaat  (iher  flmctioo  it  used  such 
diat  no  Leonard  terms  appear.)  to  adtftion,  die  average  nor¬ 
mal  velocity  at  the  rarfrice  it  now  dnqily  lelaied  to  the  total 
flow  rate  and  area,  whkh  reduces  die  kwal  jet-like  valuet  of 
waU  normal  velocity  by  dm  open  area  ratfo  ^  the  toifaoe. 
while  the  average  sUa  frioioa  it  dw  weighted  sum  of  the 
toUd  turtocc  pito  open  turtoee  eon^bodont.  Flaaffly.  the 
doubly  averaged  itieamwisr  velocity  and  turbulenoe  at  dw 
wall  are  now  nonacro. 

It  appears  dM  dw  tocadOB  la  dw  bomdary  layer  at  eddeh 
<d*d*>*0  (srhare  d  it  toe  toaeavetagii  aonaal  vetocky) 
deltoee  a  level  that  divides  toe  niad^  weakly  sheered 
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outer  flow,  which  is  readily  charaaehzable  from  even  linear 
theory,*  from  an  inner  flow  where  the  <«*!;*>  decay  results 
from  the  discrete  jet  mixing  and  streamwise  defleciion. 
Given  some  relationship  between  <«•(,•*>,.  <«•»;'>„ 

and  the  wall  structure  (where  '  denotes  turbulent  fluctuations 
and  w  wall  conditions),  the  inner  layer  typified  by  high 
shear,  relatively  low  speed  flow  may  also  be  characterizabie.’ 

In  a  constant  pressure  integral  formulation  (which  may  be 
most  convenient  for  data  analysis),  the  aforementioned 
phenomena  may  be  expressed: 

d9  /  <u>.\  (i;>.  (iTV).  C, 

dx  °  V  u,  J  u,  ~  ul  ~  ui  *2 

where  9  is  the  momentum  thickness,  Cf  the  skin  friction 
coefficient,  u,  the  freestream  velocity,  and  the  additional 
terms  correspond  to  the  three  effects  discussed  here.  («•*>,» 
having  been  deleted.  The  modeling  problem  that  the  data  in 
Ref.  I  should  significantly  impact  is  twofold: 

1)  <«>,.  and  all  depend  on  the  details 

of  the  flow  through  the  particular  wall.  Although  <ii>,  and 
<u'v">,  may  be  estimated  from  the  data,  must  be 

treated  parametrically  until  a  compatible  solution  is 
obtained. 

2)  As  in  all  TBL  modeling,  Cy  depends  explicitly  on  all 
other  effects  on  the  right-hand  side  of  the  equation. 

Therefore,  some  effort  is  required  to  assess  which  effects 
are  important  and  their  quantitative  contribution  in  a  general 
case. 

Due  to  the  relatively  low  blowing  rates  used  in  Ref.  1 
(relative  to  that  which  would  cause  separation),  both  the  vec¬ 
toring  effect  <ii>,  and  the  explicit  dependence  of  Cy  on  the 
mass  transfer  may  be  ignored  in  this  particular  case.  The 
resulting  error  of  a  few  percent  is  well  within  the  precision  of 
the  data.  This  allows  us  to  focus  on  the  representation  of  the 
terms 

-  <u’0’),/u;-{u'v">,/ul  (3) 

which  are  herein  combined  into  the  following  single  model¬ 
ing  expression 

C,  fCjb*  (4) 


rig.  I  CoMpartoM  beiwtca  the  experlweatal  data  af  Rtf.  1  and  the 
aMdtitag  of  Eq.  (3)  by  Eq.  (4). 


where 

bw  <*;>,/«,  (5) 

By  systematically  varying  the  constants  C^  and  C.,  the 
data  may  be  fit  reasonably  well  with  C,  =0.002  and 
C;  =  -  80.0  for  the  titanium  wall,  whereas  C,  =  -  0.0022  and 
C.  =  -50.0  are  indicated  for  the  sintered  wall,  as  shown  in 
Fig.  1.  Clearly,  C,  scales  with  Cf  and  C.  scales  with  C," 
implying  that  the  friction  velocity  is  the  appropriate 
reference.  Note  that,  since  the  coefficients  are  negative,  it  is 
evident  that  the  effect  of  -<ii*ti’>,  is  dominant  since  we 
expect  to  be  positive. 

If  this  functional  dependence  persists,  analysis  suggests 
that  a  rise  in  the  measured  Cf  is  possible,  over  a  range  of 
larger  b  values  than  used  in  Ref.  I,  prior  to  a  falloff  at  still 
higher  mass  transfer  rates,  it  also  should  be  mentioned  that, 
in  measurements  of  forces  on  such  a  porous  plate,  the  drag 
would  reflect  only  the  solid  surface  portion  of  the  C,  in  this 
form  of  the  boundary-layer  equations.  However,  the 
measured  force  would  include  a  pressure  difference  within 
the  passages  of  the  porous  wall  due  to  the  vectoring 

(seen  in  the  boundary-layer  equation),  which  also  involves 
the  contribution  to  Cy  averaging  from  the  open  surface. 

For  cases  in  which  the  induced  mass  transfer  vectoring, 
and  (u>.,  is  significant,  angle  (a)/magnitude  /mj  variables 
may  be  appropriate.  Utilizing  the  following  definitions  of  a. 
m: 

(6) 

m«(l -^a*)"<u>./u,  (7) 

the  vectored  blowing  term  in  the  boundary-layer  equation 
becomes 

[m* -mfl +a*)  •  1  ( -O' (I -^a*)]  (8) 

The  magnitude  m  responds  to  the  pressure  difference  across 
the  porous  surface  and  its  "loss  coefficient,"  whereas  a  pro¬ 
vides  a  sensitive  measure  of  the  induced  vectoring,  which 
also  depends  on  the  flow  within  the  porous  wall  through 
and  (a'l-'),.  The  mass  transfer  term  in  this  form 
exhibits  extrema  with  which  the  effect  of  vectoring  may  be 
evaluated.  For  a  given  a,  this  term  exhibits  an  extremum 
when  = '/i  and,  for  a  given  m,  when  <«>,/«,= 

(I  ±(1 +8ni*)'^)/4»  — m*,  Since  <ii>,  will  be  >0 

yet  <Vi,  the  effect  of  naturally  induced  vectoring  could  be 
relatively  small.  The  effect  of  actively  imposed  optimal  mass 
transfer  vectoring  is  currently  being  studied  for  both  steady 
and  unsteady  flows  and  will  be  the  subjea  of  a  future  paper. 
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Boundary  Conditions  for  Flow  Over  PenneaMe 
Surfaces 


G.  J.  Hokcason' 


Introduction 

The  subject  of  this  work  is  viscous  flow  adjacent  to  a 
permeable  surface  and  the  boundary  conditions  appropriate 
thereto  [1].  Our  objective  is  to  generalize  and  extend  the 
analysis  in  references  [2  and  3]  from  first  principles  in  order  to 
clarify  the  relationship  between  the  experimental  reality  of  a 
distribution  of  discrete  nonuniform  jets  in  a  crossflow  and 
computations  which  employ  continuous  (locally  uniform) 
normal  velocity  and  zero  slip  boundary  conditions.  The 
hypothesis  is  that  some  details  of  the  (vortical)  flow  passing 
through  the  boundary  and  the  external  (partially  slipping  (1]) 
flow  along  the  boundtuy  affect  the  qualitative  struaurc  of  the 
solution. 

The  analysis  presented  here  is  an  extension  of  the  work 
described  in  reference  [4]  and  is  directed  at  continuum  flow 
phenomena.  Although,  for  most  problems  of  practical  im¬ 
portance,  we  are  generally  interested  in  turbulent  flow,  a 
laminar  flow  is  considered  here  to  expose  the  particulars  of 
the  (»emise.  Some  of  the  resulu,  therefore,  bear  directly  on 
current  interest  in  laminar  flow  control  by  suction.  Tlie  role 
which  turbulent  flow  considerations  play  shall  be  discussed 
subsequently.  Finally,  the  restriction  that  a  boundary  layer 
structure  be  maintained,  forcing  limitations  on  the  mau 
transfer  magnitude,  shali  not  be  invoked. 


Consider  first  the  detailed  picture  of  a  porous  wall  at  a 
distribution  of  discrete  jets.  (The  cate  of  bknring  win  be  used 
for  purposes  of  discussion.)  By  placing  a  “jump-type*' 
conutri  volume  at  the  exit  plane  of  an  individual  jet,  the 
resulu  of  reference  |4)  indicate  that  the  extenial  flow  “slips** 
over  the  jet  Q-c-.  the  matt  transfer  vcaor  tiht)  with  a  vdo^ 
determined^: 

d(eiiv)  wdr,  (1) 

where  6  signiTiet  the  jump  in  values  of  the  variaMes  across  the 
coDUol  volume.  Conventional  notation  assigns  p,  u,  v,  and  r 
to  the  fluid  density,  streamwise  velocity,  normal  vdocity  and 
shear  streu,  respectivdy.  Invoking  continuity: 

»(pv)-0,  Q) 
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across  the  control  volume  allows  us  to  write  equation  (1)  as: 

5tt=6T/<w.  (3) 

This  result  shall  now  be  generalized  to  an  arbitrarily  per¬ 
meable  surface  for  comparison  to  the  aforementioned  zero 
slip,  (locally)  constant  normal  velocity  boundary  condition 
theoretical  f^ormulation.  For  cases  in  which  u  and  r  within  the 
boundary  are  zero,  equation  (3)  is  compatible  with  the  results 
in  (1-3). 

Dcvelopncal 

Consider  the  flow  along  a  porous  surface  composed  of  a 
distribution  of  openings  which  is  specified  statistically  in 
terms  of  size  and  location.  The  extern^  flow  is  assumed  to  be 
two-dimensional,  however,  extension  of  the  results  to  include 
three-dimensional  effects  is  straightforward  and  will  be  noted. 
A  local  averaging  process  is  now  invoked  [3]  in  planes  parallel 
to  the  boundary.’  The  spatial  scales  of  the  averaging  are  large 
relative  to  the  porosity  yet  small  relative  to  the  streamwise 
(and  span-wise)  variations  in  the  mean  flow.  Each  of  the 
variables  is  decomposed  into  its  local  spatial  mean  plus  a 
perturbation  (e.g.,  u  a  <u>  +  u").  When  this  decomposition 
is  applied  to  the  Navier-Stokes  equations  for  constant  density 
flow,  the  resulting  equations  may  be  spatially-averaged  to 
provide  the  result: 

<«>,  +  <t;>,=0  (4) 

(<u><u>),= 

+  (5) 

<t<>,  +  (<«i><tO),-t-(<i;>=),= 

-<P)  ,/P-<X>, -<M“t'“>r -*’<«>>  (;u.  (*) 

where  subscripts  denote  differentiation  and  subscripts  in 
parentheses  denote  a  vector  component.’  Not  surprisingly, 
this  is  reminiscent  of  Reynolds  decomposition  for  turbulent 
flow.  Additional  terms  required  for  three-dimensional  mean 
flows  are  apparent.  If  we  now  apply  the  jump-type  control 
volume  to  the  permeable  boundary  over  the  averaging 
domain,  the  results  analogous  to  reference  {4|  are: 

5<f>=0  (7) 

<i;>6<u>= -5<u''t;®>  +  i-5<«)(j,  (8) 

/O 

5(5/)^)=- (9) 

In  this  representation,  total  mass  flow  through  the  per¬ 
meable  boundary  and  shear  force  on  it  are  locally  equivalent 
(over  the  averaging  area)  between  the  discrete  and  continuous 
representations.  The  average  mass  flux  is  reduced  from  that  in 
the  discrete  jets  by  the  ratio  of  open  to  total  surface  area.  The 
average  shear  stress  is  reduced  from  that  on  the  solid  in¬ 
terstices  by  the  ratio  of  solid  surface  area  to  total  surface  area. 
Some  new  features  appear,  however.  Most  importantly  the 
average  streamwise  velocity  at  the  boundary  is  nonzero  and 
related  to  that  in  the  boundary  (allowing  for  the  possibility  of 
obliquely-oriented  jets)  by  equation  (8).  The  subtleties  of  such 
a  situation  (whether  due  to  geometrical  or  flow-induced 
obliquity)  regarding  effects  of  internal  normal  stresses  versus 
surface  shear  stresses  and  interpretation  of  experiment-theory 
comparisons  with  the  appropriate  boundary  conditions  is 


*  Note  that  ihc  averasinf  weightinf  function  was  chosen  10  be  a  comtani  over 
ihe  averagini  domain  and  lero  outside  of  it.  In  this  manner  Reynolds  averaiing 
rules  are  retained  and  the  tolled  Leonard  terms  do  not  appear. 

’This  formulation  also  signifieamly  clarifies  the  effect  of  man  transfer  on 
the  wall  vcrtiaiy  flux.  If  <a>  is  nonzero  fdue  10  active  vectoring  or  naturally 
with  nomsnally  normal  injection),  <ia>i;i,  responds  to  the  four  advective terms 
in  enuaihm  (S>  even  in  consiam  pressure  flows.  If  <■>  is  also  nonsleady,  the 
<«),  term  also  modifies  the  vorticity  flux  at  ihc  wall. 


clarined  by  equation  (8).  This  relationship  also  embodies 
information  on  the  injea^  flow  vorticity  and  structure. 

Applying  the  same  decomposition,  averaging  and  jump- 
type  boundary  control  volume  analysis  to  the  full  equations  in 
vorticity  variables  results  in: 

+  i>8<ci.'>,;,,,  and  (10) 

(11) 

which  provides  often-needed  information  on  flow  curvature 
at  the  wall. 

Discttssioo 

Therefore,  we  have  established  a  relationship  between  the 
discrete  distribution  of  jets  in  experiments  and  the  continuous 
mass  flux/shear  stress  distributions  in  computational 
predictions,  as  long  as  <u>#0  at  ihe  boundary.  If  (u)  is  set 
equal  to  zero  at  the  boundary,  calibration  of  some  aspeas  of 
the  theoretical  formulation  with  experimental  data  must 
compensate  for  the  approximation.  It  may  be  that  the 
magnitude  of  the  slip  is  small,  inherently  or  due  to  term 
cancellation  in  equation  (8).  However,  the  structure  of  the 
injected  flow  does  now  enter  the  problem.  Possibly,  the 
appropriate  slip  velocity  (guided  by  equation  (1))  may  be 
obtained  empirically  by  judicious  iteration  between  ob¬ 
servation  and  theory. 

However,  the  formulation  suggests  several  more  systematic 
approaches.  First,  since  jet  structures  may  be  reasonably 
parameterized,  the  situation  is  much  simpler  than  turbulence 
in  representing  the  nonlinear  terms.  Therefore,  a  quantitative 
assessment  of  the  effect  of  differences  in  detailed  structure 
between  porous  boundaries  which  are  equivalent  in  the  mean 
is  possible.  Secondly,  it  suggests  that  applying  computational 
boundary  conditions  at  some  height  off  the  wall  for  which 
<u®v®>  a>  0,  analogous  to  the  flow  over  a  wall  film,  may  be 
more  appropriate. 

Consider  these  two  uses  of  equation  (8).  Across  the 
boundary,  5<u>  may  be  computed  if  we  know  something  of 
the  injected  flow  structure.  It  is  clear  that,  over  the  openings 
u"  and  t’®  >  0  and  over  the  solid  portion  u®  and  f®  <  0. 
Therefore,  (m®i;®>  is  a  positive  number.  Given  <v>,  <u>  at  the 
wall  may  be  computed  (as  a  function  of  <ui>  at  the  wall,  as 
required  computationally)  by  specifying  typical  detailed 
distributions  of  u®  and  u®  of  interest  within  the  porous 
boundary,  computing  their  jump  across  the  boundvy  ac¬ 
cording  to  equation  (1)  and  then  averaging.  The  procedure  is 
necessarily  iterative,  requiring  an  initial  guess  of  the  slip 
velocity. 

Furthermore,  assuming  that  the  height  above  the  boundary 
at  which  (u®!?®)  >  0  is  small,  the  jump  expressions  may  be 
applied  (to  some  lesser  degree  of  accuracy)  between  it  and  the 
waJi  boundary.  Note  that  the  jump  across  the  boundary  must 
first  be  computed.  In  this  manner  the  equations  are  used  to 
provide  the  boundary  conditions  «ti>  and  <u>)  for  com¬ 
putations  at  that  level  where  the  discrete  jets  have  mixed  out. 
This  is  reminiscent  of  using  wall  funaions  and,  indeed,  some 
profile  hypothesis  is  required  to  do  this  accurately,  (f  explicit 
account  of  the  injected  turbulence  is  of  interest,  the  tur¬ 
bulence  kinetic  energy  and  length  scale  at  the  displaced  height 
could  be  computed  from  model  equations  utilizing  the  jump- 
type  control  volume  analysis.  This  is  relevant  to  the  comment 
in  reference  [3]  in  which  the  authors  refer  to  blowing  as 
"provoking  the  viscous  sub-layer."  Utilizing  the  approach 
developed  here,  the  extent  to  which  the  injected  flow  structure 
totally  disrupts  the  near  wall  conditions  may  be  assessed. 

Consider  the  two-equation  kt  model  as  an  example  only, 
not  necessarily  propos^  here  as  a  general  computational  tod 
for  these  flows.  In  this  preliminary  evaluation,  we  assume 
that  the  low  Reynolds  number  postulates,  relative  to  the  near 
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wall  region  of  the  flow  o\er  an  impermeable  boundary  {6.  ■?]. 
are  inappropriate  here  and  the  relevant  form  of  the  equations 
IS  that  for  which  molecular  effects  are  ignorable.  By  applying 
the  jump  conditions  across  the  layer,  6v,  in  wrhich  the  jets 
"mix  out.”  the  following  expressions  for  the  turbulence 
kinetic  energy  and  isotropic  dissipation  ensue: 

=  6^. -r  IP, -Di  |6y,  and  (12) 

pi6f  =6q  +  IP  -ZJ,  lov,  (13) 

where  I  I  indicate?  the  vertical  spatial  average  across  the  5> 
layer  and  the  diffusive  fluxes  (q),  production  (P)  and 
dissipation  {D)  are: 

q.  =(p./a,  )dk,'d\  ,  q,  )dt/dy 

P.  =p  (dii/di  )*'  ,  P,/(  =  C^  ic  P^/k  (14) 

D,  =<  ,  D,/t  =  CixDi,/k 

Utilizing  expressions  or  values  for  n,  (C^pk^/t),  a,.  C,, 

and  C; .  the  magnitudes  of  k  and  t  at  by  m.ay  be  computed  in 
terms  of  k  and  <  injected  through  the  boundary  and  the  near 
wall  profiles.  Therefore,  the  formalism  by  which  the  external 
flow  responds  to  the  strucutre  of  the  injected  flow,  which  was 
observed  earlier  for  nonturbuleni  vortical  details,  carries  over 
to  turbulence  as  well 

Note  that,  for  flows  with  sireamwise  pressure  gradients, 
injection  of  material  which  undergoes  a  phase  change  at  the 
boundary  complicates  the  discussion  of  injected  vortical 
structure  due  to  the  V<p>  x  V<p>  '  vorticity  generation 
effect  In  addition,  dilatation  effects  may  be  important  if 
either  phase  is  compressible.  Finally,  for  integral  for¬ 


mulation,  the  standard  mass  transfer  term  is  modified  to:  (I 
-  <M>„/u,)  (<!>./«,)  and  (u®!’”)  and  the  Reynolds 
stresses  (if  the  flow  is  turbulent)  at  the  wall  are  included. 
Utilizing  such  a  formulation,  the  effects  of  slip/vecioring  and 
the  peculiarities  of  various  walls  on  <u®i;®>  may  be  examined 
with  data  such  as  that  in  reference  |8]. 
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Vortidty  with  Variable  Viscosity 
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The  formulations  developed  in  Ref.  I  are  eminently 
uKful  in  the  analysis  of  voniciiy  mechanics  in  comptea 
flows.  As  noted  therein,  the  focus  of  the  study  was  not 
viscous  compressible  flow.  However,  some  of  the  equations 
main  the  viscosity  inside  various  derivatives,  impiying  that 
an  the  variable  visGoaity  effects  are  accounted  for.  at  least 
formally.  The  puipose  of  this  Note  is  to  document  some  ef¬ 
fects  of  the  additional  terms  required  to  complete  the  for¬ 
mulation  when  variable  viscosity  flow  is  the  subject  of 
interest. 

Following  the  notation  in  Ref.  2.  the  equations  of  motion 
may  be  written  in  the  following  form:  4 

where  e„  it  the  symmetric  pan  of  the  deformation 

(lu./Air.  *  ka./kx, ) 

T — ‘ - 


Awr^  and  6^  it  the  Kronccker  delta.  Density,  vonicity, 
velocity,  normal  sircu,  body  force,  and  the  coefTidcni  of 
viscosity  are  assigned  their  traditional  aomeaclatatt. 


(I) 


(2) 
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Straightforward  differentiation  of  the  viscous  terms  in  E<). 
(1)  allows  them  to  be  written 


"r  j  ^ - 3 - - 


(3) 


In  addition,  the  appropriate  vector  calculus  identity  allows 
this  expression  to  be  rewritten  as 


-m(  V  > 


4  iidA  /  \ 


(*) 


In  order  to  derive  the  vorticity  transport  equation,  this  ex¬ 
pression  is  substituted  into  the  equations  of  motion,  each 
term  is  divided  by  density  p,  and  the  entire  expression  is  sub¬ 
jected  to  the  curl  operator.  Comparing  this  with  Eq.  (I)  in 
Ref.  I.  it  is  clear  that  the  additional  variable  viscosity  terms 


P 

must  be  addressed  for  general  viscous  flows. 

The  focus  of  interest  in  this  study  is  the  effect  of  the  first 
term  in  flows  for  which  the  approximation  pa  constant, 
AaO  is  a  good  one,  yet  n  must  be  allowed  to  vary.  In  addi¬ 
tion  to  various  laminar  flows,  this  situation  is  applicable  to 
turbulence  modeling  as  well,  in  which  p  is  related  to  the  local 
flowfield  properties.  We  consider,  therefore,  the  following 
term  in  the  general  vorticity  transport  balance: 

For  cases  of  interest  in  which  «i  varies  predominantly  in  one 
direction,  herein  Xj.  the  components  of  the  curl  in  the  three 
orthogonal  directions  are 


^  1 

ix,  ’ 

k  iXy  ) 

'  ix,  ' 

x^ 

^  ( 

1  e, yin) 

f  filing. 

If. 

ix,  ' 

^  iXy  / 

ix, 

<  ix,  r 

•*1 

^  I 

feyyin) 

f  ei;3>tY 

Jf, 

ix,  ' 

^  iXi  / 

axj 

^  ix,  )' 

If  X| ,  Xi.  and  Xy  are  associated  with  the  streamwise,  nomul, 
and  spanwise  coordinates  of  a  shear  flow,  it  is  clear  that  ex¬ 
plicit  three-dimensionality  is  required  to  impact  streamwise 
voniciiy  from  variable  viscosity  terms.  However,  in  two- 
dimensional  flow,  the  spanwise  vorticity  is  affected  by  the 
component  of  the  curl  in  x,: 


dxy  V  dx,  dxi  ) 

which  may  be  written 


(8) 


(9) 


If,  in  addition.  Uy^u,,  this  expression  further  reduces  to 
one-half  of 

ill  /  d  iu,  \  i^ii  in, 

ix,  \3xj  ixj  /  djrJ  3xj  (10) 


in  i^u,  i*ii  iu, 
dxj  ixi  ixi 


(II) 


If.  for  example,  n  exhibits  some  empirical  dependence  on, 
say.  iu,/ixj,  the  vonicity  generation  will  respond  to  a  term 
proponionaJ  to 


which,  clearly,  can  undergo  several  sign  changes  across  a 
shear  flow  and  contribute  positive  or  negative  terms  to  the 
vorticity  balance.  The  accurate  representation  of  variable 
viscosity  effects  is,  therefore,  critical  in  representing  the 
detailed  vortical  mechanics. 
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Jet  Nozzies 
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Recently  interest  in  thrust  veaor  control  with  rigidly 
mounted  thrusters  has  been  revived.  The  results  described 
in  this  Note  are  associated  with  the  side-force  generated  by  a 
supersonic  jet  with  an  exit  plane  which  is  oblique  to  the 
centerbne  of  the  exhaust  flow.  As  shown  in  Fig.  I,  depending 
on  whether  the  flow  is  over-  or  under -expanded,  the  direction 
of  the  exhaust  flow  is  deflected  either  toward  or  away  from  the 
normal  to  the  exit  plane  by  an  amount  related  to  the  pressure 
imbalance  and  the  projeaed  area  of  the  exit  plane  in  the 
lateral  direction.'  Therefore,  control  over  the  direction  and 
amount  of  side-force  is  available  by  varying  the  stagnation 
pressure  and  exit  plane  obliquity. 

Although  the  present  results  are  associated  with  a  two- 
dimensional  conTiguration,  with  the  exit  plane  obliquity 
csublished  by  extending  one  wall  of  the  nozzle,  application  to 
an  axisymmetric  configuration  is  conceivable.  For  example, 
consider  a  fixed  axisymmetric  thruster  enclosed  by  a  cylin¬ 
drical  housing  which  is  split  longitudinally  into  various  sec¬ 
tions  (see  Fig.  1).  At  the  null  position,  the  end  of  the  cylinder 
coincides  with  the  exit  plane  of  the  thruster.  As  various  com¬ 
binations  of  cylindrical  sections  are  translated  forward  of  the 
exit  plane,  the  radially  projected  area  of  the  exit  plane  may  be 
oriented  toward  any  azimuthal  angle.  The  amount  of  pro¬ 
jeaed  area  and  its  orienution  may  also  be  changed  rapidly, 
which  adds  another  control  feature  to  the  thrust  veaoring. 

The  focus  of  the  work  described  in  Ref.  2  involved  the  ex¬ 
tent  to  which  viscous  effeas  affeaed  the  performance,  par¬ 
ticularly  at  over-expanded  conditions.  A  two-dimensionai 
supersonic  ja  of  aspea  ratio  (height/ width)  equal  to  1.5  was 
established  over  the  Mach  number  range  1 .0  to  3.0.  One  wall 
of  the  nozzle  contour  was  then  extended  a  specified  fraaion  of 
the  nozzle  height  to  establish  an  oblique  exit  plane. 
Transparent  and  parallel  lateral  sidewalls  bounded  the  nozzle 
and  extensions  thaeof  through  which  the  flow  angle  could  be 
observed.  The  side  force  generated  was  measured  with  a  strain 
gage  balance  to  an  accuracy  of  ±2^  of  the  reading,  based  on 
the  scatter  in  repeated  experimenu  and  comparison  to  the 
measured  flow  angle. 

In  order  to  assess  the  degree  to  which  viscous  effecu  impede 
accurate  analytical  prediaion,  the  amount  of  asyninwtry  was 
Umited  to  maintain  the  wall  extension  within  the  dtHnain 
ciMblished  by  the  Mach  wave  oiiginatiag  from  the  opposite 
nozzle  wall  at  the  exit  plane. 

Figures  2-4  depia  the  forces  measured  on  three  nozzles 
whose  exit  plane  is  offsa  by  an  amount  1.0,  4/3,  and  5/3 
timet  the  jet  height,  respectively.  Comparison  with  predictiont 
of  two-dimensional  inviscid  flow  theory  is  shown  also.  These 
resuhs  may  be  converted  into  force  cocfflcicnu  by  dividing  the 
rcMilu  by  the  surface  area  of  the  wall  extension  and  the  exit 
dynamic  pressure.  The  corresponding  plot  of  the  Fig.  3  dau  it 
presented  in  Fig.  5.  These  remits  indiwe  that,  even  in  the  sim¬ 
ple  ease  estabttshed  here,  oontriviai  deviationt  between  in- 
visdd  flow  theory  and  experiment  arc  obaerved,  although  the 
general  trends  predicted  by  analysis  are  coma.  It  is  also  dear 
that  the  veaoring  control  may  extend  well  imo  the  over- 
enpanded  region  before  the  nozzle  boundary  layers  begin  to 
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sepwttc  and  the  entire  flow  degradei.  ahhoiigh  tlw  Reynold* 
number  dependence  of  the  effect  in  this  conflgumtion  remaio* 
#  to  be  (luantified.  On  the  baiis  of  tbeae  data,  it  is  condudsd 

that  control  of  the  thruster  ptenun  pressure  and  exit  plane 
obliquity  (including  the  aforementiooad  axisymmetric  case 
wherein  the  skteforce  vector  may  he  oriented  azimuthally)  pro¬ 
vides  the  ability  to  quickly  orient  the  thrust  vector  of  nongim* 
balled  jets. 
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DR  APT  OF  FUTURE  PUBLICATION 


The  angular  structure  of  viscous  incompressible  flow 

Part  I.  Analysis 

Gustave  J.  Hokenson,  Chief  Scientist 
The  Hokenson  Company 
Los  Angeles,  CA  90005 

Abstract 

The  governing  equations  for  time-dependent  three-dimensional  viscous  shear  flows  are 
transformed  into  velocity  vector  angle  and  magnitude  variables  to  expose  the  flow  angle 
structure  in  various  limiting  situations.  For  purposes  of  analysis  when  the  mean  stream- 
wise  flow  is  reasonably  well  characterizable,  this  formulation  admits  the  possibility  for 
extracting  essential  flowfield  peculiarities  and  providing  insight  into  various  complex 
interactions,  including  those  associated  with  turbulent  flow.  For  full  numerical  simulations 
of  shear  flows,  the  forrixilation  simplifies  the  mathematical  structure  and  may  allow  for 
more  effective  computation. 


I.  Introduction 


The  use  of  velocity  vector  magnitude  and  angle  variables  in  the  solution  of  inviscid 

1  2 

compressible  and  incompressible  flow's  has  a  long  history  of  successful  application.  ’ 

Less  common  is  the  implementation  of  such  an  approach  in  viscous  shear  flows,  which  is 
the  subject  of  this  study.  Part  I  of  the  work  involves  an  exposition  of  the  general  formu¬ 
lation  and  its  application  in  various  limiting  cases.  In  a  future  paper,  the  results  of  a 
computational  study,  which  has  been  initiated,  shall  be  presented. 

Consider  the  time-dependent  three-dimensional  continuity  and  Navier-Stokes  equations 
for  a  constant  density  Newtonian  fluid^: 


+  V  +  w 

y  2 

=  0 

9 

(1) 

Ut  - 

+  (uv)y  ^ 

(uw)^ 

=  F 

"t  " 

(vu)^ 

>  (vMy  * 

( vw)^ 

=  G 

(2) 

Wt  . 

(wu)x 

+  ( wv )  + 

y 

(w*)^ 

=  H 

f 

where: 


F  =  -  p^  ♦  vV*u 

G  =  -  Py  v7*v  (3) 

H  =  -  ♦  v7*w  , 


-B3- 


and  V  is  the  constant  kinematic  viscosity,  the  conventional  Laplacian  operator  and  p 
is  defined  as  the  mean  normal  stress  (herein  equated  to  the  equilibrium  thermodynamic 
pressure)  divided  by  density.  Following  conventional  terminology  u,x;  v,y;  w,2  correspond 
to  the  nominal  streamwise,  normal  (along  which  the  dominant  velocity  gradients  occur), 
and  spanwise  velocities  and  coordinates,  respectively. 

The  following  dependent  variables  are  now  introduced  in  an  attempt  to  expose  the  angu¬ 
lar  structure  of  various  shear  flows: 

m  E  (u’  +  V*  w*  )/u^  ^ 

e  =  v/u  (4) 

a  =  w/u  , 

where  is  an  appropriate  reference  velocity.  In  general,  we  do  not  expect  this  trans¬ 
formation  to  be  useful  when  the  streamline  lies  in  the  (y-z)  plane  orthogonal  to  the 
nominal  streamwise  (x)  direction.  That  is,  when  u=0  but  v,w  ^  0  we  may  have  to  retreat 
from  this  formulation,  locally. 

Clearly,  m  may  be  expressed  as: 

m  =  (u/u^j,j)*  (1  +  0^  ♦  a*)  .  (5) 

For  convenience  in  transforming  the  equations,  several  of  these  terms  are  grouped 
together  such  that  u  may  be  written: 


-BA- 


1 


u^=mB;6"‘5(l  +  e^+a^)/u^j.^  .  (6) 

Utilizing  tnese  definitions  and  relationships  for  m,  0  and  a,  Eqs.  (3)  become: 

(mB)^^  +  (nn6)^  +  (mse)^  +  (mBa)^  =  F  (7) 

((mB)^9)^  +  (mSe)^  +  (mgeMy  +  (mBaO)^  =  G  (8) 

((mS)^a)j.  +  (mBa)^  +  (mBae)y  +  (mSaM^  =  H  .  (9) 

By  expanding  the  differentiation  in  Eqs.  8  and  9  and  utilizing  Eq.  7,  the  following  com¬ 

pact,  yet  complete,  relationships  for  G  and  a  ensue: 

(10) 

(11) 

where: 

Njj,  =  3/3t  -  u  O/ax  +  09/9y  +  a9/9z)  .  (12) 

Note  that  Eqs.  10  and  11  may  be  combined  to  express  the  sum  and  difference  of  the  two 
flow  angles  in  the  following  form: 

(0  t  a)  =  |(GiH)-  F(0  ±  a)}/u  .  (13) 
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Cleariy,  continuity  and  the  magnitude,  m,  momentum  equation  must  be  considered  in 
order  to  provide  a  closed  set  of  equations  for  u,v,w  and  p.  They  shall,  indeed,  be  add¬ 
ressed  later  in  this  analysis  but  the  focus  of  the  v^ork  in  Part  I  is  to  expose  the  structure 
of  a  and  o  for  various  cases  in  which  u  can  be  characterized  parametrically.  This  approach 
embodies  the  dynamics  of  interest  in  lieu  of  the  simple  kinematics  derivable  from  contin¬ 
uity  if  sufficient  information  on  the  streamwise  and  cross-stream  derivatives  is  given. 

Prior  to  addressing  various  specific  cases,  the  form  of  the  governing  equations  encourages 

us  to  examine  the  loci  where  u  =  0.  As  long  as  the  flow  remains  incompressible,  this 

mathematical  structure  is  complete  and  G,  a  and  their  derivatives  are  everywhere  bounded 

within  the  previously-stated  restriction  of  streamlines  from  the  y-z  plane.  (It  has  been 

found  that  several  important  cases  of  interest  satisfy  this  restriction  when  analyzed  in 

detail  but  do  not  in  traditional  formulations.  An  example  is  a  boundary  layer  with  wall 

mass  transfer.  Clearly,  at  the  wall  v  0  but  u  is  traditionally  set  equal  to  zero.  It  is  now 

4  5 

known  that  u  at  the  wall  is  small  but  finite’.) Therefore,  when  u  =  0,  =  9/9t,  0  =  G/F 

and  a  =  H/F  with  the  values  of  (G  -  0F)/u  and  (H  -  aF)/u  determined  from  L 'Hospital's 
rule.  If,  in  addition,  u  =  u^  =  0,  as  occurs  at  the  point  of  separation  in  a  wall-bounded 
shear  flow,  G  -  9F  and  H  -  aF  values  are  stationary  at  that  point  (still  possibly  zero) 
and  the  values  of  (G  -  0F)/u  and  (H  -  aF)/u  are  determined  from  the  curvatures  of  G  -  0F 
and  H  -  oF.  For  various  cases  G  -  0F  and  H  -  oF  may  be  either  maxima  or  minima  and, 
therefore,  3a/9t  and  90/St  may  be  positive  or  negative.  Subsequently,  we  shall  analyze 
specific  steady  flows  which  modify  this  argument  somewhat. 

II.  Discussion 


We  now  consider  some  restricted  cases  of  the  governing  equations  in  more  detail. 


-B6- 


A.  Steady  Twc  -  Dimensional  Flow 


(14) 

or,  in  streamfunction  (\|!,  x)  variables: 

=  (G  -  eF)/u*  .  (15) 

Within  the  same  constraints  as  before  (0  »  by  restriction  and  0,0  »  by  assumption 

X  y 

of  physicality),  \A/hen  u  =  0,  0  =  G/F  but  also,  G  -  0F  is  stationary  at  zero.  If,  in  addition, 

Uy  =  0,  G  -  0F  and  its  first  three  derivatives  must  be  zero.  This  is  clearly  a  significantly 
strengthened  constraint  from  its  unsteady  flow  counterpart. 

Utilizing  the  definitions  of  F  and  G  along  with  the  angle/magnitude  transformation,  the 
0  equation  may  be  written: 

(1  -  2v(mB)J/m8)0^  +  (0  -  2v(m8)  J/mB)©„  =  (-p„  +  ©p^)/mB  +  vV^OArr^)^.  ( 16) 

In  order  to  isolate  certain  cases  of  interest,  it  is  convenient  to  non-dimensionalize  this 
equation  with  length  scales  and  L^yj,  angle  scale  0^^.^  and  the  previous  velocity  scale 
u^j.y  utilizing  the  following  definitions: 


In  tnis  situation,  a  =  H  =  0  and  the  governing  0  equation  becomes: 

0^  ^  GGy  =  (G  -  0F)/u* 
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X  = 

‘"(y)'''-(x) 

III 

C 

c 

M 

(  X ) 

X 

®(y) 

III 

N) 

C 

c 

M 

y 

Re  E 

''(r)'-(x) 

e  =  Re‘* 


the  governing  0  equation  in  dimensionless  variables  becomes: 

"  -  “  ‘®Px  ’  Py''^®  ( r )  * 

'<®xx  *  ®yy^^n/u  .  (18) 

In  order  to  retain  the  two-dimensional  equation  of  continuity  in  its  primitive  form,  X  is 
set  equal  to  0^^.^  so  that  the  0  equation  may  be  written: 

(1  -  .  (6  -e»(y)/xM0y  =  (ep,  -  Py/X’)/u>  *  e(0„  *  0yy/ 

A“)/u  .  (19) 


we  now  consider  various  cases  of  interest  for  which  e  «  1  .  The  most  straightforward 
limit  is  X  »  1  for  which  the  first  order  solution  for  ©  (in  an  implied  perturbation  expan¬ 
sion;  may  be  developed  by  letting  X  -►  ®  in  Eq.  19,  resulting  in: 

0^  +  00y  =  0Pj^/u*  ♦  e  ,  (20) 

which  can  be  made  formally  invariant  with  respect  to  e  by  utilizing  x/c  and  y/e  as 
Independent  variables. 
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If,  on  the  other  hand,  X  «  1  such  that  e/A^  «  1,  i.e.  X  is  small  relative  to  one  but  stiil 
large  relative  to  the  presumed  ve;y  small  e  ,  the  governing  0  equation  becomes: 


0,  *  eey.= 


•eyy/U- 


(21) 


where: 


u '  =  Xu  ,  and 

£’  =  e/X  , 


(22) 


are  used  to  maintain  the  canonical  form  of  Eq.  20  which  is  independent  of  X.  Evidently, 
the  appropriate  velocity  in  this  case  is  a  typical  v.  Eq.  21  also  may  be  made  invariant 
with  respect  to  t  via  x/e'  and  y/e*. 


This  case  is  to  be  contrasted  with  the  traditional  boundary  layer  limit  for  which  X  «  1 

but  e/X^  =  0(1)  (i.e.  L  /L  ~  Re“^)  so  that  all  the  coefficients  in  Eq.  19  are  of  order  unity. 

y  ^ 

If  e  =  0(x*) ,  the  term  -  0^  must  be  retained  and  the  equation  becomes: 

0^  +  (0  -  (c/X^)  ♦(y))ey  r  -  Py/X^u*  +  (e/xM0yy/u  .(23) 

For  purposes  of  analysis,  is  treated  as  a  small  (variable  viscosity  -  like)  term  which 
can  be  incorporated  as  a  perturbation  correction  to  the  solution  of: 

+  00y  =  -  Py/U'*  ♦  C*0yy/U'  ,  (24) 
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which  is  essentially  identical  to  Eq.  21,  except  that  c'  =  e/X  =  .  Eq.  2A  also  may  be 

made  invariant  with  respect  to  e  by  employing  x/e^  and  y/c^  as  independent  variables. 

The  traditional  use  of  Eqs.  21  and  24  involves  the  formal  limit  X  -►  0  for  which  p^  0, 
the  boundary  layer  assumption.  Note  that  the  behavior  in  this  limit  is  formally  the  same 
as  our  previous  analysis  of  u  =  0. 


Our  use  of  the  previous  scaling  is  directed  at  simplifying  the  coefficients  in  the  governing 
equation  while  retaining  all  the  terms  such  that  the  mathematical  structure  is  complete. 
We  leave  it  to  the  solution  to  expose  the  behavior  of  the  dependent  variables  when  c  and 
X  are  small  enough  to  have  allowed  various  terms  in  the  equation  to  be  neglected. 


Further  analysis  of  Eq.  24  is  facilitated  by  the  use  of  streamfunction  independent  variables 
(ii»,x)  wnich  allows  the  equation  to  be  written: 


Sx  =  - 


(25) 


If  we  define  an  equilbrium,  0^^^  which  is  independent  of  x,  its  distribution  in  terms 
of  u'  and  p  is: 


e^0(eq.)  =  /l/u'/dp/u*  dij;  +  ♦  C2 

which,  if  u'  varies  weakly  with  ^  (as  it  must  in  order  to  account  for  ♦ 
theory),  reduces  Eq.  (26)  to: 


(y) 


(26) 

via  perturbation 
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^^®(eq.)  “  (Jpctv^)/u’*  +  +  C2 


(27) 


In  addition,  thin  internal  layers  may  exist  within  the  shear  flow  for  which  Eq.  (25)  provi¬ 
des  the  following  'jump'  (6)  condition: 


Jdp/u'  =  e^6(u'0^) 


(28) 


or,  for  weakly  varying  u': 


g’6G^  =  6p/u ' * 


(29) 


Finally,  the  general  solution  of  Eq.  (25),  in  terms  of  a  known  or  characterizable  p  and 
u'  field.  Is  readily  obtained  from  transient  heat  conduction  solutions  to  the  equation 
written  as: 


Gy  -  e^(u'e^)^  =  -  p^/u' 


(30) 


and  the  u'  (vl;)  variations  treated  perturbatively  and  parametrically.  In  addition,  since  the 
u'  (\l;)  variations  are  weak  in  order  to  treat  via  perturbations,  u'  may  be  set  equal 
to  a  constant  determined  by  collocation  in  order  to  optimize  the  approximation  of  neg¬ 
lecting  its  distribution. 

If  u'  (0)  characterizations  are  meaningful,  Eq.  30  may  be  written: 


where: 


S'x  -  ®'w  '  ■  %  ’ 


u'd0  =  d©' 


(31) 
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which,  although  non-linear,  is  particularly  convenient  for  analysis.  The  'equilibrium'  solu¬ 
tion  ana  'jump'  conditions  of  Eq.  (31)  are  apparent. 

Analysis  of  the  behavior  of  Eq.  20  is  more  clearly  carried  out  computationally  inasmuch 
as  the  conventional  streamf unction  (\l,x)  transformation  results  in  a  complex  elliptic 
equation  for  0.  If  the  less  conventional  (\l',y)  streamfunction  transformation  is  invoked 
0/5x  =  -  v3/?t;  o/5y  =  3/9y  +  uo/Sijj),  Eq.  20  may  be  written  in  terms  of  0*/2  =  ^  : 

%  '  -  .  (52) 

in  which  the  non-linearities  limit  closed  form  analysis.  However,  the  equilibrium  and 
'jump'  conditions  may  be  approached  as  before,  particularly  if  u  (^ )  relationships  (where 
ud?  =  d0')  are  useful. 

B.  Steady  Three  -  Dimensional  Flow 

By  expanding  and  non-dimensionalizing  Eqs.  10  and  11,  the  following  expressions  for  ©  and 
a  ensue  as  c  0: 


-B12- 


where  the  additional  scale  ratio  introduced  and  the  characte¬ 

ristic  ancle,  has  been  set  equal  to  in  order  to  retain  continuity  in  its  primitive 
form. 


In  the  lifTiit  as  X^^^  -►  ®  ,  Eqs.  (33)  may  be  written: 

{?/dx  >  G5/8y  +  a3/3z  -  P^^/u*  -  (e/u)9V9xM  e*a  =  0  .(34) 

Alternatively,  if  X^^^  =  X^^^  -*>  0  such  that  e/X‘  is  still  small,  Eqs.  33  become: 

{9/9x  +  e9/9y  +  a9/9z  +  -  ( e  ' /u '  )  ( 9  V  9y  *  +  9V92*)}0,a  =  0, 

(35) 

where  j  =  y  for  G  and  j  =  z  for  a  .  In  the  traditional  boundary  layer  limit  this  same  form 
is  retained  with  replacing  e'  and  the  (presumed  weak)  variable  viscosity  -  like  $ 
terms  accounted  for  perturbatively.  Clearly,  if  and  are  large,  they  must  be 
included  with  0  and  a  in  the  advective  operator. 

In  the  mixed  limit  case,  X^^^  -►  <*>  and  \y) "  °  such  that  e/X*^y^ 

'^3®  =  -  Py/u'*  +  ^yy'"^’ 

Njtt  =  OPj^/U*  +  ®yy''*J 

where:  =  9/3x  +  09/9y  ♦  a3/92  , 


=  1  ,  Eqs.  (33)  become: 


(36) 


-B13- 


and  must  be  accounted  for  at  higher  order.  Alternatively,  if  -*■  ®  and  -*•  0 
such  that  e/X’.  =  1  ,  the  Eqs.  (33)  are  v/ritten: 


N^e  = 

N^a  =  -  p^/u'*  +  e^a^^/u' 


(37) 


again  assuming  that  is  subsequently  accomodated.  By  scrutinizing  these  mixed  limit 
cases  it  is  clear  that  the  essential  angular  structure  of  the  flo\«  will  be  markedly  diffe¬ 
rent  from  that  we  have  observed  previously.  Consider,  for  example,  the  case  where  p^  = 
0  and  a  is  small.  This  results  in: 


®x  "  ®®y  =  ®Z2^^  ’ 

indeed  a  curious  mix  of  derivatives.  Three-dimensional  phenomena  are  most  clearly  inve¬ 
stigated  numerically  and  will  be  described  In  a  future  publication.  Note  that,  unlike  two- 
dimensional  flows,  total  knowledge  of  u  is  inadequate  to  prescribe  the  angular  structure 
from  the  kinematics  of  the  incompressible  continuity  equation. 


C.  Time  -  Dependent  Flow 


By  adding  the  dimensional  terms: 


0^/u  and  a^/u  (39) 

to  the  advective  operators  in  Eqs.  (33),  the  three-dimensional  formulation  embraces 
unsteady  flows  as  well.  After  non-dimensionalization,  these  terms  are  scaled  by  a 
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coefficient; 


^(x)^*-(r)''(r  ) 


(40) 


wnere  reolaces  and  a^ter  clearing  the  equations  with  respect  to 
and  respectively.  If  t^^.^  is  chosen  as  additional  dimensionless 

terms  ©^./u  and  a^/u  are  present  in  all  equations  for  which  X  -*•  «»  . 

In  cases  for  which  -*■  0  ,  it  is  nxire  consistent,  formally,  to  set  t^j.^  = 

and  express  the  operators  as: 


a^/u'  +0^^+ . 


,  and 


(41) 


For  mixed  limits,  the  a  and  0  equations  respond  most  naturally  with  different  'time  con¬ 
stants',  significantly  enriching  the  range  of  solution  characteristics  to  be  exposed  numeri¬ 
cally. 

For  the  two-dimensional  flow  governed  by  Eq.  (24),  the  unsteady  counterpart  is: 


0t/u'  0^  ♦  GOy  =  -  Py/u'*  ♦  €*0yy/u' 


I '  *  +  €^( 


(42) 


or,  in  conventional  streamfunction  variables: 


Gj  u'0^  =  -  P^;  ♦  (u'G^)^ 


(43) 
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If  trie  transformation  =  u'  dO  allows  for  meaningful  parameterization  of  the  velocity 
field,  Eq.  (43)  may  be  written: 


(44) 


both  of  which  illustrate  significant  changes  in  the  boundary  layer  approximation  p^  =  0. 


For  purooses  of  analysis,  if  u'  =  u'  (li/),  along  curves  in  the  x  -  t  plane  for  which  x  -  u't 
is  a  constant,  Eq.  43  may  be  written: 


0t  =  -  p^  +  e^u*  (u'G^)^  .  (45) 

Once  again,  the  general  transient  structure  is  so  complex  for  either  two-  or  three-dimen¬ 
sional  flows  as  to  require  numerical  investigation. 

with  regard  to  time-dependent  flows,  considerable  interest  exists  relative  to  the  angular 
structure  of  turbulent  flows.  The  counterpart  of  the  Reynolds  stress  tensor  may  be  found 
by  setting  e  =  0  in  Eqs.  (33).  For  notational  simplicity,  the  case  =  X^^j  =  *  will  be 
used  without  justification.  In  this  case,  the  equations  for  0  and  a  become: 

{u‘*a/8t  +  a/3x  +  03/3y  +  ot3/3z  -  a,©  =  0  .(46) 

Since  the  continuity  equation  may  be  written; 
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Eqs.  (46)  become: 


a.  *  (ua)  (u9a)  +  (uaM  =  ap^/u  ,  and 

L  y  A 

G^,  +  (uG)^  (uOMy  (uae)^  =  ep^/u 


(48) 


If  we  consider  two-dimensional  "turbulence"  in  order  to  expose  the  angular  Reynolds 
stresses,  Eq.  48  may  be  written  utilizing  traditional  Reynolds  decomposition  as: 


®x  ®®y  ^  -  0(P’u’  +  (u'0'  )y )/u  -  IPT’‘(20y  +  P)/u  -  0^'  , 

(49) 

where  P  =  Pj^  and  triple  correlations  of  the  fluctuations  (designated  by  the  prime)  and 
the  X  -  derivative  of  have  been  neglected.  In  addition,  the  ensemble  mean  values 
(designated  by  the  overbar)  have  been  assumed  to  be  time  -  independent.  The  term  (i7F')y 
is  reminiscent  of  the  conventional  Reynolds  shear  stress.  However,  this  formulation  expo¬ 
ses  additional  explicit  effects  of  P,  0  and  and  correlations  between  P',  u*  and  ©'.  It 
remains  to  be  seen  whether  or  not  the  explicit  appearance  of  these  various  terms  is 
sufficiently  advantageous  in  light  of  the  additional  superficial  complication. 

With  regard  to  the  turbulent  burst  dynamics,  various  aspects  of  the  vortical  instabilities 
are  often  described  in  terms  of  a  four  quadrant  u’  -  v’  map.  ^  Equations  for  0'  and  its 
statistical  characterization  may  also  be  derived  from  the  full  0  and  0  relationships,  provi¬ 
ding  an  analytical  tool  with  which  to  quantify  the  various  details  which  support  the  afore¬ 
mentioned  four  quadrant  uniquenesses  of  turbulent  burst  dynamics. 
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III.  Conclusion 


The  full  time-dependent  three-dimensional  Navier  -  Stokes  equations  for  incompressible 
viscous  shear  flows  have  been  converted  to  vector  angle  and  magnitude  variables  for 
analysis  in  various  limiting  cases.  This  dependent  variable  transformation,  in  conjunction 
with  a  streamfunctioh  independent  variable  is  seen  to  simplify  the  system  of  equations 
and  permit  analysis  of  the  angular  structure  with  the  pressure  and  velocity  field  treated 
parametrically.  The  utility  of  this  formulation  in  the  study  of  turbulence  appears  to  have 
merit.  Finally,  the  equation  structure  is  often  considerably  more  "symmetric"  in  these 
variables  and  may  facilitate  numerical  studies  which  are  currently  under  way. 
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TBL  INSTABILITY  COMPUTATION  INPUTS  AND  RESULTS 


Figs.  Cl  &  C2.  Input  Mean  Velocity  and  Eddy  Viscosity  Profiles  and 

Their  Derivatives. 
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Fig.  C3  &  Table  C3.  Alpha  Solutions  For  The  Particular  Case  Shown. 
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Fiy.  C3 


SOLUTIONS  OF  ALPHA 


BETA 
OMEGA 
EPSMAG 
TAUEPS 
RE  NO. 


3.00 


0.100 

2.000 

1.000 

0.000 

44369 


2.00 


THE  HOKEMSON  COHPANY 

Table  C3 


BETA  e.iee 

ONEGA  2.000 
EPSNAG  1.000 
TAUEPS  0.000 
RE  NO.  44369.000 


ALPHA  REAL 

ALPHA  INAGINARY 

K 

PHASE  SPEED 

-0 . 30270 

-1.37840 

4.55368 

-6.60720 

0.13009 

-1 .06779 

-8.2081 1 

15.37404 

-0 . 25528 

-0.79761 

3.12447 

-7.83455 

-0.78500 

- 1 . 65000 

2.35669 

-2.54777 

0 . 82936 

-0.71168 

-0.8S811 

2.41150 

1.06350 

-0.58185 

-0.54711 

1.88058 

1 .28000 

-0 . 46000 

-0.35936 

1 .56250 

1.42000 

-0 . 36000 

-0.25352 

1.40845 

1 .54000 

-0.34000 

-0.22078 

1.29870 

0.52391 

-0.89511 

-1.70851 

3.81744 

-1.19534 

-2.49390 

2.01106 

-1.67316 

-1.74600 

-2.70400 

1.54868 

-1.14546 

-0.09086 

9.64849 

-7.13706 

-22.01140 

-1.36535 

1 .53441 

-1.12382 

-1.46463 

Figs.  C4  -  CIO.  Eigenfunctions  For  The  Solution:  ot^.  =  .13  and  =  -  1*07  On  Fig.  C3. 
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Fig.  C9 
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CIO 

TYZ  STRESS  PERTUBATION 

REAL  IMAGINARY 


Figs.  &  Tables  Cll  -  C13.  Alpha  Solutions  For  The  Particular  Cases  Shown 
With  Conditions  Identical  to  Fig.  C3  Except  is  Non-Zero. 
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Table  Cll 


I 


BETA  0.100 
OnEGA  2.000 
EPSHAG  1.000 
TAUEPS  1.000 
RE  NO.  44369.000 


ALPHA  REAL 

ALPHA  I HAG I NARY 

K 

0.04822 

-0.61343 

-12.72114 

-0.73597 

-1 .69037 

2.29679 

-1.16810 

-2. 12207 

1 .77120 

0.38616 

-1 . 16745 

-3.02318 

0.00584 

-1.38015 

-236.44127 

-0.44789 

-1 .51448 

3.38136 

0 . 68070 

-0 . 96358 

-1.41556 

0.68239 

-0.44285 

-0.64897 

0.90013 

-0.81632 

-0.90690 

1.11133 

-0.71862 

-0.64663 

1.18127 

-0.00043 

-0.00036 

1 .27817 

-0.59046 

-0.46196 

1 . 38630 

-0.50984 

-0.36777 

-0 . 34523 

0 . 30375 

-0.87984 

-0.48100 

0.70200 

-1.45946 

-0.95000 

2.20000 

-2.31579 

I 


» 


PHASE  SPEED 

41 .47530 
-2.71749 
-1 .66932 
5. 17613 
342.63199 
-4.46537 
2.93813 
2.93086 
2.22190 
1.79965 
1 .69310 
1 .56473 
1 .44268 
-5.79327 
-4.15800 
-2.10526 
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Fig.  C12 
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Table  Ci2 


BETA 
OnEGA 
EPSnAG 
•  TAUEPS 

RE  NO. 


0. 100 
2.000 
1  .000 
10.000 
44369.000 


ALPHA  REAL 

ALPHA  I HAG I NARY 

K 

PHASE  SPEED 

0.10000 

-0 . 46000 

-4.60000 

20.00000 

-0.18000 

-1.32000 

7.33333 

-11.11111 

0.10000 

-1.20000 

-12.00000 

20.00000 

0.58000 

-1 . 12000 

-1.93103 

3.44828 

0.90000 

-0.93000 

-1 .03333 

2.22222 

1 . 14000 

-0.78000 

-0.68421 

1 .75439 

1.32000 

-0 . 68000 

-0.51515 

1.51515 

1.52000 

-0.60000 

-0.39474 

1 .31579 

1 . 17000 

-0 . 20000 

-0.17094 

1 .70940 

-1 .07429 

2.35657 

-2.19360 

-1 .86169 

-0.02175 

0.00801 

-0.36826 

-91.94980 
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Fig.  C13 


SOLUTIONS  OF  ALPHA 


BETA  0.100 
OMEGA  2.000 
EPSMAG  1.000 
TAUEPS  30.000 
RE  NO.  44369. 


3.00 


2.00 


1.00 


0.00 


-1.00 


0.00 


1.00 


2.00 


ALPHA  REAL 
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Table  Cl  3 


BETA  0.100 
OnEGA  2.000 
EPSRAG  1.000 
TAUEPS  30.000 
RE  NO.  44369.000 


ALPHA  REAL 

ALPHA  I HAG I H ARY 

K 

PHASE  SPEED 

0.06000 

-0.44000 

-7.33333 

33.33334 

-0.64000 

- 1 . 72000 

2.68750 

-3. 12500 

-1.26000 

-2.16000 

1 .71429 

-1 .58730 

0.60000 

-1 . 10000 

-1 .83333 

3.33333 

1.00000 

-0.90000 

-0.90000 

2.00000 

1.20000 

-0 . 76000 

-0.63333 

1 .66667 

1 .20000 

-0.20000 

-0.16667 

1 .66667 

1.36000 

-0.66000 

-0.48529 

1 . 47059 

0.14000 

-1 . 18000 

-8.42857 

14.28571 

-0.12000 

-1.30000 

10.83333 

-16.66667 

-1 .11325 

2.40235 

-2.15796 

-1.79654 

-0.00350 

-0.00050 

0.14286 

-571.42853 

C22 


Figs.  Cl^  a  &  b  and  Table  C14.  A  Modified  Input  Eddy  Viscosity  Profile  With  The 
Peak  Located  Twice  As  Far  From  The  wall  As  In  Previous  Solutions.  Alpha 
Solutions  For  The  Particular  Case  Shown. 
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Table  Cl^ 


BETA  0.100 
OnEGA  2.000 
EPSnAG  1.000 
TAUEPS  1.000 
RE  NO.  44369.000 


ALPHA  REAL 

ALPHA  IHAGINARY 

K 

PHASE  SPEED 

-0.98400 

-2.02950 

2.06250 

-2.03252 

-0.46000 

-1 .66000 

3.45833 

-4.16667 

-0.19600 

-1 .33600 

6.74747 

-10.10101 

-0.19800 

-0.54600 

2.75758 

-10.10101 

0.17400 

-1.38000 

-7.93103 

1 1 . 49425 

0.51500 

-1 . 14427 

-2.22187 

3.88350 

0 . 75600 

-1 .01800 

-1 .34656 

2.64550 

0.98400 

-0.91450 

-0.92937 

2.03252 

1.16250 

-0.80500 

-0.69247 

1 .72043 

1 . 34830 

-0.70360 

-0.52184 

1 . 48335 

0.88800 

-0.41640 

-0 . 46892 

2.25225 

1 . 23800 

-0 . 40960 

-0 . 33086 

1.61551 

1.41900 

-0.16600 

-0.11698 

1 .40944 

1 . 18700 

-0.00400 

-0.00337 

1 .68492 

-0.31800 

-0.03200 

0.10063 

-6.28931 

-0.20234 

0.59344 

-2.93289 

-9.88435 

-0 . 98600 

2.26225 

-2.29437 

-2.02840 

